Undoing measurement-induced dephasing in circuit QED 
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We analyze the backaction of homodyne detection and photodetection on superconducting qubits 
in circuit quantum electrodynamics. Although both measurement schemes give rise to backaction in 
the form of stochastic phase rotations, which leads to dephasing, we show that this can be perfectly 
undone provided that the measurement signal is fully accounted for. This result improves upon 
that of Ref. [1], showing that the method suggested can be made to realize a perfect two-qubit 
parity measurement. We propose a benchmarking experiment on a single qubit to demonstrate the 
method using homodyne detection. By analyzing the limited measurement efficiency of the detector 
and bandwidth of the amplifier, we show that the parameter values necessary to see the effect are 
within the limits of existing technology. 
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I. INTRODUCTION 

Recent developments in the field of superconducting 
qubits realizing long qubit coherence times [2], non- 
destructive single-shot homodyne measurements [3], use 
of fast FPGA's [4, -^j], and approaching quantum lim- 
ited parametric amplifiers [()], have opened up for the 
realization of e.g. entanglement generation by measure- 
ment [7, S] and the implementation of error correction 
protocols and fault tolerant quantum computing [()]. It 
is known that measuring a quantum system introduces 
backaction which sometimes can limit the fidelity of the 
measurement. To realize a high-fidelity measurement, it 
is therefore necessary to understand the nature of the un- 
necessary backaction and how it can be mitigated or, in 
the best case, avoided. 

In this paper, we discuss how to completely undo the 
dephasing induced by measurement on qubits in a circuit 
quantum electrodynamics (cQED) setup [10, 11]. Our 
analysis is focused on the read-out of a single qubit, but 
it is equally valid for the case of joint measurement on two 
qubits, which enables perfect parity measurement in this 
system. The reason for largely limiting the discussion to 
the singlc-qubit case is that the effect can be equally well 
understood and observed in a one-qubit measurement, 
which is beneficial from an experimental point of view. 

In circuit QED, the qubit is coupled to a microwave 
resonator which gives rise to a state dependent shift of 
the resonator frequency. The qubit and the field states 
become entangled and it is possible to read out the state 
of the qubit by measuring the state of the emitted mi- 
crowave field. The homodyne measurement of one and 
two qubits was analyzed in Refs. [12, 1.3] and [14], respec- 
tively. In the latter it was shown that when measurement 
parameters are tuned to realize a parity measurement, 
there is unwanted backaction in the form of a stochastic 
phase causing dephasing in the post measurement state. 
In Ref. [I] it was shown that this backaction could be 
partly undone by recording the homodyne current and 
applying conditional control pulses to undo the phase 



uncertainty. 

In this work, we extend the analysis of Ref. [1] to show 
that not only part of the measurement-induced dephas- 
ing, but all of it, can be undone provided that the mea- 
surement signal is turned off and all photons pumped 
into the resonator arc accounted for. The result is in 
agreement with the simpler picture of a strong projective 
measurement of one of the field quadratures for which it 
is easy to see that the dephasing can be undone given that 
no knowledge about the initial superposition is obtained. 

To further explore the connection between information 
gain and dephasing, we consider the backaction when the 
field is measured using direct photodetection. This is in- 
teresting from an experimental point of view due to the 
recent developments in photon counting devices in the 
microwave spectrum [1-5-17]. Theorywise, the read-out 
is interesting because the measurement adds additional 
backaction apart from the fundamental projection pos- 
tulated by quantum mechanics [In]. By deriving an ef- 
fective stochastic master equation for the qubit degrees 
of freedom, we show that, when the photons contain no 
information about the qubit state, it is possible to com- 
pletely undo the dephasing by recording the time when 
each photon is detected. Hence, the connection between 
information and backaction is established in this case as 
well. 

Finally, we suggest the implementation of a simple one- 
qubit experiment to test the undoing of dephasing in the 
homodyne case and show that it is possible to observe the 
effect with existing state of the art measurement devices. 

This paper is organized as follows. In Sec. II we show 
the connection between information gain and backaction 
using a simple toy model consisting of two entangled 
qubits and a projective measurement. In Sec. Ill we 
present the model describing measurement of one qubit in 
a resonator, considering two cases: homodyne detection 
and photodetection. In Sec. IV we show, for both detec- 
tion schemes, how the measurement-induced dephasing 
can be perfectly undone. For homodyne measurement, 
the result holds equally well for the case of two qubits. 
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In Sec. V we describe an experiment to test this in cir- 
cuit QED for the case of homodyne detection, and discuss 
suitable parameters. Conclusions are given in Sec. VI. 



II. INFORMATION GAIN VS. BACKACTION 

To illustrate the connection between information gain 
and backaction we consider two entangled qubits, 

\i^)=a\gg)+l3\ee). (1) 

It is clear that making a measurement along the z-axis 
of qubit two, <jf'\ gives information about |ap and |/3p, 
with the post measurement state of the first qubit given 
by the mixed state pqbi = diag(|Q;p,|/3p). If we, on the 

(2) 

other hand, choose to measure cri , we will get the results 
±1 with equal probability, hence revealing no information 
about a and j3. The post measurement state of qubit one 
is 

_{a\g)+P\e), if ai'' = +1, 
\ "Iff) - /3|e), if cri ' -1, 

resulting in a pure state, which is either left untouched or 
phase flipped. Note that, depending on the measurement 
result, one can reverse the effect of measurement, leaving 
the state uncorrupted. 

The situation in cQED is very similar to this model 
albeit more complex in two senses; the second qubit 
is replaced with the more complex field state, and the 
measurement cannot be considered strong but is instead 
happening on a timescale comparable to the resonator 
linewidth. In the following we show that the simple sym- 
metry between information gain and backaction is pre- 
served despite these complications. 



III. SYSTEM 

A. Hamiltonian and master equation 

We consider a circuit QED setup with one qubit cou- 
pled to a resonator driven by a measurement signal de- 
scribed by the Hamiltonian [10] 

H = ~^'^z + i^rO^ CL + g {a) + a(T_|_) 

+ {ed{t)a^e-''^^' + e*,{t)ae''^^'), (2) 

where we have set h = 1, ujaiyJr) is the resonance fre- 
quency of the qubit (resonator) , g is the coupling strength 
between the resonator and the qubit, a{a)) is the anni- 
hilation (creation) operator for the resonator, (t_((t+) is 
the lowering (raising) operator of the qubit, £d{t) is the 
amplitude of the measurement drive, and lo^ is the fre- 
quency of the measurement signal. 



In the dispersive regime, |A| = jwa — Wr-I ^ g, the 
effective Hamiltonian of the system is given by [10] 

^off= '^(J^ + /\ra^a + xa^aG^ + {£d{t)a^ + £*d{t)a) , (3) 

where x — 5^/ A, >^a = "^a+Xi — w^— w^, and we have 
moved to a frame rotating with the measurement signal 
frequency tud- The master equation for the cavity-f-qubit 
system in this frame is [!!)] 

P = -i [Hcs, p] + kV [a] p + "fiV [(T_] p 

+ (4) 

where p is the density matrix of the system, 2? [X] p = 
X pX'^ - ^X'^ X p- ^pX^ X arc Lindblad operators [20], k 
is the resonator decay rate and 71 , 7^ are the qubit relax- 
ation and dephasing rates respectively. If the resonator 
state is coherent, and k 3> 71, the field amplitudes solv- 
ing Eq. (4) are given by [21] 

"9(0 = -i£dit) -i{Ar- x) oig{t) - ■^ag(<), (5) 

de(t) = ~i£d{t) - i (A,. + x) ae{t) - ■^ae(i), (6) 

depending on whether the qubit is in the ground state 
\g) or the excited state |e). 



B. Stochastic master equations 

Although the master equation in Eq. (4) describes the 
evolution of the open system, we need to take into ac- 
count the conditional evolution given one certain mea- 
surement trace. In this section, we consider two types of 
measurement; homodyne detection and photodetection. 

The stochastic master equation (SME) [22] describing 
the system evolution conditioned on homodyne detection 
is given by 

dp = CtotpAt + [ae-"!'] pAW{t), (7) 

where £totP is the RHS of Eq. (4), 77 is the measurement 
efficiency and M.[c]p = cp -\- pc^ — + c') p is the mea- 
surement operator with (A) = tr (Ap). The phase of the 
local oscillator is given by cj), and dW{t) is a Wiener in- 
crement defined by E [dW{t)] = and E [dW^(i)2] = dt. 
The homodyne current is given by 

j{t)dt = y/lWT{ae-"t' + a^e'-^) dt + dW{t). (8) 

For the case of photodetection the SME is [23] 

dp = -1 [H,s, p] d< -I- (1 - r,)KV [a] pdt 
n 

+jiV[a^]pdt+^V[cT,]pdt 

+g [a] pdN{t) - ^7]kM [a^'a] pdt, (9) 
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where G [c] p — j§f^^P, and dN{t) is a Poisson point pro- 
cess with the properties dN{t)^ = dN{t) and E [dN{t)] = 
(a^a) dt. 

Although Eq. (7) and Eq. (9) can be used directly to 
study the measurement backaction, they provide little 
intuition and understanding. Instead we choose to trace 
out the resonator degrees of freedom and work with ef- 
fective equations for the qubit degrees of freedom only. 



C. Effective stochastic master equations 

The effective SME for the qubit degrees of freedom 
when read out through homodyne detection was derived 
in Ref. [12] and is given by 



dp 



,qb 



B 



2 

70 + r 



dt + ^iV[a-]p^''dt 



-V [a,] p'i^dt 



UqM [n„e-"^] p'^''dW{t), 



(10) 



where IIq. = ctglig + aeRe- Here Td = 2xlm (aga*) is the 
measurement-induced dephasing and B = 2xRe(aga*) 
gives the AC Stark shift. The homodyne current is given 
in terms of qubit operators by 

j{t)dt ^ ^{iV^e'^^ + n^e^*)^,, dt + dW{t). (11) 

In the Appendix, we follow the procedure in Ref. [12] to 
derive an effective SME for the case of photodetection. 
Applying the displacement transformation 



ngD[ag]+n,D[ae]. (12) 



where 11, 



\&/g){e-/g\ and D [a] = exp {aa^ — a' 



a) 



Eq. (9) takes the system to a frame of reference where 
the resonator field is in the vacuum state. This makes 
it possible to trace out the resonator degrees of freedom, 
enabling us to derive the effective SME for the qubit con- 
ditioned on photodetection 



dp' 



B 



2 

70 + Td 



[(T„p<ib] dt + ^^V[a^]p'^''dt 
P [cr^] /9*dt - [Ha] p^'^'dt (13) 



^i^kM [nln„] p^i^dt + g [n„] p^^'^dNit). 



IV. DEPHASING 

To study the unwanted measurement-induced dephas- 
ing in the single qubit case, we consider the case = 0, 
which implies ag = —a* for the field state. In this case, 
doing homodyne detection with (j) = tt/2 (measuring the 
projection of ai on the imaginary axis) or photodetection 
(measuring the photon number) will not give any infor- 
mation about the qubit state. However, the measurement 



still gives rise to dephasing in the form of stochastic phase 
rotations. For the homodyne case, the analysis is equally 
valid for a measurement of two-qubit parity where the 
unwanted backaction gives rise to dephasing within the 
negative parity subspace []]. In analogy with Sec. II, we 
show how this can be perfectly undone given the mea- 
surement result. 



A. Homodyne detection 

We consider the effective SME, Eq. (10), with 71 = 
70 = 0, and exclude the coherent evolution since this only 
gives rise to a deterministic phase which can be undone. 
In this case the SME can be solved analytically with the 
solution 



Pggit) = PggiO), 

Peg{t) = Peg (0) exp 



~Td{s) + 2KriRc^ {ag{s))) ds 



+2iy/J^ J Re{ag{s))dWis) ), 



(14) 



where, from now on, we drop the superscript "qb". The 
second equation gives the dephasing and was investigated 
in Ref. [].] for the two-qubit case in the context of par- 
ity measurement. There, a constant measurement tone 
was assumed and it was found that there was inevitable 
dephasing resulting in a mixed post-measurement state. 
This was attributed to the finite response time of the 
cavity. Here, we extend the analysis and turn off the 
measurement signal at time imoas such that ed{t) = 
e?ri (0(0 ~ ~ ^mcas)), whcrc Q{t) is the Heaviside 
step function. In this case the solution to Eq. (5) is given 
by 



Ugit) 



2e, 



2x + in 
1 — exp 



1 — exp 



ix)t))e{t) (15) 



(y- - iXj {t - imoas)) ^ Q{t - tmcas) 



which inserted in the real part of the exponent in Eq. (14) 
gives 



lim* 



dim* 



'Td{s) + 2KRe^ {ag{s))) ds 



iag{t) - a,{t)f 



0, 



(16) 



where it should be noted that we have assumed a quan- 
tum limited measurement i] ~ 1. Hence, by turning off 
the measurement signal and allowing for all the informa- 
tion that is encoded about the qubit state to be mea- 
sured, we completely recover the initial pure state of the 
qubit except for the stochastic phase seen in the second 
row of Eq. (14). This is one of the main results of this 
paper. If, on the other hand, the measurement signal 
isn't turned off, there will always be some photons left in 
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the resonator and the information which is thus withheld 
from us gives the residual dephasing that was found in 
Ref. [1]. 

The fact that the measurement induced dephasing is 
zero when not acquiring information about the qubit 
state is not apparent when looking at Eq. (10). However, 
since the decrease of the off diagonal element is associated 
with the gradual projection onto one of the measurement 
eigenstates, the result is expected. This was also found 
in [13] using a quantum bayesian formalism in the limit 
X < K. 

The imaginary part of the exponent in Eq. (14) gives 
stochastic phase kicks to the qubit. Averaging this over 
many trajectories would give rise to dephasing. However, 
since we measure the homodyne current, which simplifies 
from Eq. (11) to 



j{t)dt = 2y/]^m{ag{t)) dt + dW{t) 



(17) 



for the case (j) — 7''/2, and we know Q:g{t), we can extract 
dW{t) and the stochastic phase kicks from our measure- 
ment. In conclusion, if = 1 we are able to completely 
determine the measurement-induced phase kick to the 
qubit and undo it after each run of the experiment, essen- 
tially undoing the measurement-induced dephasing. The 
elimination of the stochastic phase was also discussed in 
[13]. Here, a direct feedback controller was used to mod- 
ulate the qubit frequency. This suggestion requires the 
modulation to be done on the time-scale of the fluctu- 
ations in the measured current. Our suggestion is not 
limited by this since all the phase is undone using a sin- 
gle control pulse at the end of the measurement. 

The connection between phase rotation and measure- 
ment result can be made intuitively clear by consider- 
ing the strong version of the measurement, i.e. a pro- 
jection on one of the eigenstates of the field quadrature 



P = i(at 



qubit is given by 



i) I The post measurement state of the 



|V')qb = {6{p\ag)\g) + -f{p\ae)\e)) , 



(18) 



where TV is a normalization and the overlap integral 
is given by [24] 



1 



ri/4 



exp 



-Im(/3)2 - 2/32 - (p + ^^/2/3)^ 



and 5 and 7 are arbitrary complex numbers satisfying 
|(5|^ + ItI^ — 1 and P\p) = p\p). The off-diagonal element 
in the post measurement density matrix is hence given 
by 



post 



-i2\/2Rc(Qg)p 



(20) 



where we have used the specific symmetries of the mea- 
surement setup. In the strong measurement limit, the 
only effect on the state is to apply a stochastic phase ro- 
tation between \g) and |e) whose magnitude depends on 
the measurement outcome p. This is in analogy with the 



analysis above. The difference between the strong and 
weak measurement only lies in the fact that the entire 
measurement trace must be accounted for if the mea- 
surement is weak. 

In Fig. 1(a) we plot the purity as a function of the wait- 
ing time after the measurement signal has been turned 
off, toff. Using the procedure outlined above to correct 
the measurement-induced phase kicks, we see that the 
post measurement state approaches a pure state when 
all the photons have been accounted for. Observe that 
the purity oscillates during the measurement such that 
the qubit state is pure at certain points in time. We can 
understand this by looking at the time evolution of the 
resonator states which is plotted in Fig. 1(b). After the 
measurement signal has been turned off ag and ae spi- 
ral towards the origin. We see that the times for which 
P = 1 coincides with the times at which the resonator 
states overlap. The overlap integral in Eq. (19) is the 
same at these points such that the photons remaining in 
the resonator carry no information about the qubit state. 
To ensure that the purity remains 1, however, we must 
wait until all photons have left the resonator, which takes 
a few resonator lifetimes. 



Parity measurement 

The result presented in Eq. (16) can easily be extended 
to the case of joint parity measurement on two qubits 
which was studied in detail in Ref. [1]. In this case, 
the qubit-resonator couplings and detunings are given by 
91 = 92 and Ai = — A2 which enables an observer to 
distinguish the parity of a two-qubit state. As shown in 
[14], this results in dephasing within the negative par- 
ity subspace spanned by \eg) and |(?e). In analogy with 
the single-qubit measurement, this dephasing originates 
from the entanglement with the field states a^e and a^g . 
Since these evolve in the same way as ag and ae (but 
with the replacement x — > 2^ in Eqs. (5) and (6)), the 
result in Eq. (16) holds equally well with the substitu- 



tion ag — >■ age and ae — > aeg, which shows that the 



measurement-induced dephasing can be completely un- 
done in the two-qubit case as well. In this way it is pos- 
sible to realize a perfect parity measurement in cQED. 



B. Photodetection 

In Sec. IV A, we saw that it is possible to gain insight 
about the backaction of the weak measurement by look- 
ing at the corresponding strong measurement limit. In 
the case of photodetection the detector adds additional 
backaction to the system by absorbing the photons that 
are detected. Hence, the post-measurement state is vac- 
uum and not an eigenstate of the photon number oper- 
ator corresponding to the number of measured photons. 
For the choice of parameters considered in Sec. IV A there 
is however nothing in the field states which allows us 
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iDh 0.9- 




ing and has the analytic solution [25] 
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Figure 1. (Color online) a) The purity as a function of the 
waiting time toa- After each measurement, a phase correc- 
tion conditioned on the measurement trace was applied as 
described in the text. The parameters are given by Ed = n, 
imcas = 5/k and x ~ The inset shows the purity as a 
function of measurement efficiency for the case of; no feed- 
back (red, dotted), feedback without accounting for all pho- 
tons (blue, dashed) and feedback with full accounting (solid 
green). Here, ton — 3/k. 

b) The time evolution of the cavity field Ofe (green) and Qg 
(red, dashed). After the measurement signal is turned off, the 
field starts to spiral back towards the origin. At certain points 
in time (arrows) = ctg, for which the qubit state is pure. 



to distinguish the qubit states when measuring photon 
number. We now show that, despite the additional back- 
action, the purity of the qubit state is preserved as it 
should be if no information is acquired. To this end we 
consider the effective photodetection SME, Eq. (13), with 
the same choice of parameters as for the homodyne case 



—iikM [n^n,] pdt + g [n,] pAN{t). (21) 

This gives the following two equations for the elements 
of the density matrix: 



<^Pgg 
dpeg 



0, 

{-Td + VK {al + Nc,)) pegdt 



N, 



No. 



-pegdN{t), 



(22) 



(23) 



where we have introduced the average photon number 
Na = \ctg\^- The second equation describes the dephas- 



Peg{t) = exp (^J^ {-Tdis) + rjK (a^(s) + Na{s))) ds 

xPes(O) n -e''«(*"', (24) 



where ^ = Arg(Q:e) and are the times at which the 
photons are detected. Observing that 

al = 2Re^ {ag) - No.+ilai{al) , (25) 

Eq. (24) reduces to (in the case 77 = 1 and t ^ tmoas) 

(0 = Peg (0) exp ( i / Klm(ag(s)) ds j -e 



Peg 



(26) 



where we have used the result from Eq. (16). We see 
that the qubit state acquires a phase kick at each time 
a photon is detected. The size of this is just twice the 
phase of the applied read out field which is a deterministic 
quantity. As for the homodyne case, the measurement 
backaction preserves the purity of the state and the initial 
state can be perfectly recovered by an observer who keeps 
track of the times that the photons arrive at the detector. 



V. EXPERIMENT 

Although experimental progress in microwave pho- 
todetectors is promising [L'^-IT], it is, as of today, difficult 
to test Eq. (26) experimentally. Homodyne detection, on 
the other hand, is routinely used to read out the single 
and multi-qubit states in cQED [2()], suggesting that, at 
least, it should be possible to undo the measurement- 
induced dephasing with this setup. The feedback scheme 
suggested above assumes a perfect record of the mea- 
surement current needed to calculate the necessary con- 
ditional phase kick back. Unfortunately, this requires 
a quantum limited {tj = 1) detector with infinite band- 
width (in order to perfectly resolve j(i)). Current state of 
the art measurements in cQED arc done with paramet- 
ric amplifiers [G, 27, 28] having a typical measurement 
efficiency of 77 = 0.4 and bandwidth BW = 10 MHz. 

To see if the effect of feedback can be observed with 
realistic measurement parameters, wc model the effect of 
finite amplifier BW by inserting a second resonator be- 
tween the system of interest and an amplifier having infi- 
nite bandwidth and detector with measurement efficiency 
T]. The second resonator is characterized by a loss rate 
Kf, = BW/2 through each of its sides and annihilation 
operator h. The setup is depicted in Fig. 2. Using the 
(S;L,H) formalism for cascaded quantum systems [29, 30] 
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Figure 2. The setup used to model finite bandwidth of the 
detector. 
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Figure 3. (Color online) The purity of the qubit density 
matrix after measurement, with phase correction after t = 
tmeas + ioff, as a fuuctiou of amplifier bandwidth. Parameters 
are as in Fig. 1 with tos — 10/k. The inset shows the purity 
vs measurement efficiency for BW — 10k. 



the SME for the full system is given by 



dp- 



[H^s, p] + kV [a] p + 2KbV [b] p 



-^,V[a_]p+^V[a,]p 



dt 



+y/J^M [he-">'] pdW{t), (27) 

where, in this case, the homodyne current is given by 

j{t)dt ^ ^/J^{be-'"^ + h^e"!') dt + dW{t). (28) 

In Fig. 3 we plot the purity of the post-measurement 
state as a function of the amplifier bandwidth. The re- 
sults are obtained by numerical integration of Eq. (27) 
including both resonator's degrees of freedom. Given 
the homodyne current, we extract the noise process by 
assuming Eq. (17). Since the true current is given by 
Eq. (28), this will lead to a wrong estimation of the phase 
kick back which limits the fidelity of the process. For the 
parameter values quoted above, there is however a clear 
distinction between the purity of the post-measurement 
state with and without the feedback protocol proposed, 
thereby suggesting that it should be possible to see this 
effect using state of the art cQED read-out and control 
technology. Such an experiment would serve as a power- 
ful benchmark to show how well qubits in a cQED archi- 
tecture can be controlled and measured. 

Apart from the measurement infidelities in the exper- 
imental setup, we have to take into account the time it 



takes to analyze the measurement result needed to calcu- 
late the right kick back. Comparing the coherence times 
of transmon type qubits (^ 50/is [2]) with current per- 
formance of state of the art FPGA's [4], this might be 
possible but comes at the cost of difficult implementa- 
tion. Instead of the direct scheme proposed in Sec. IV A, 
we therefore propose the following, simple protocol which 
avoids these complications 



1. Prepare the qubit in the superposition state 
(|0> + |1))/V2. 



2. Apply a weak "measurement" signal with local os- 
cillator phase (/) = 7r/2 during time imcas and wait 
a time ioff- Record the corresponding homodyne 
current and store the result. This will only induce 
a stochastic phase rotation which we seek to undo. 

3. Apply a pulse that changes the phase of the qubit 
randomly by 



4. 



Do a strong measurement of the qubit (0 = 7r/2). 
This does not cause additional phase rotations but 
projects the qubit state on the measurement eigen- 
basis. Store the result. 



When this protocol has been repeated many times 
there is plenty of time to calculate the required phase 
0caic for each measurement run. Given this, post-select 
the runs where 4>q is in the range [4>ca\c — S, ^icaic + <5] and 
look at the purity of the final state. Here, (5 is a phase 
interval which determines the trade-off between purity 
of the chosen states and the number of successful kick 
backs. 



VI. CONCLUSION 

We have studied the measurement-induced dephasing 
in a single qubit in the two cases of homodyne detection 
and photodetection. We show that the dephasing can be 
completely undone if all the photons in the resonator are 
measured. We also show that this procedure undoes all 
the unwanted backaction in a parity measurement on two 
qubits, thereby improving the results of Ref. [1]. Further- 
more, we have described and analyzed an experimental 
proposal for testing the effect in the case of homodyne 
detection on one qubit in a resonator. We show that the 
experiment is feasible in a parameter regime accessible 
with existing technologies in circuit QED. 
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Appendix: Deriving the effective stochastic master 
equation for photodetection 



Q with A4 [c] p = cp + pc) and Q [c\p = cpc) — p) 
df p^] At + (1 - ^)nV [a^] p^At 

+7il? [^^] p^dt + [af] p^dt 
--r]KM [a^PflP] fdt + g [flP] fdN{t) 



(A.l) 



To derive the effective SME for photodetection, 
Eq. (13), from the full SME, Eq. (9), we follow the pro- 
cedure outlined in Ref. [12] which was used to derive where we have defined operators in the displaced frame 
Eq. (10) from Eq. (7). We begin by applying the dis- = P^^cP. By a direct calculation of the transformed 
placement transformation, Eq. (12), on the unnormalized operators and using Eqs. (5) and (6) we arrive at the 
version of Eq. (9) (replacing the superoperators A4 and SME in the displaced frame 

I 



dp' 



^<j, + {Ar + x<Jz) a^a + i {edUi + e*n„) ,p^ dt + {1 - ?/) kV [a] fdt + jiV [al] fdt 
V [a,] p^dt + (1 - r;) « (2? [n„] p^ + a [p^, + [n„, p^] a)) dt 
^ijK {{a^a + 2alli + nj,n„) p^ + f {a^ a + 2a^n^ + nj,n„)) dt 
ap^a^ + ap^ni + H^fa^ + Hap^Ul - p^]dN{t). 



1± 
2 
1 



(A.2) 



r 



The transformed density matrix can be expanded in such that the reduced density matrix for the qubit, p'^''\ 



eigenstates of CTz and a^c 



is found by tracing out the resonator degrees of freedom 



trres (Pp^P^) 



(A.4) 



= E E pln.,\i){3\{m\DHa,)D{a.,)\n) 

n,m— i.j—g,e 
oo 

= E E Pnr,izj\i){j\drnnijexp{-ilm{aja*)) , 

n,m— i.'j—g^e 



where we have defined dp 
Ui — a,. 



{p\D{l3,j)\q), and ft. 



oo The evolution of the density matrix elements are given 

E E PlrmM,i){m,jl (A.3) byEq. (A.2) 



n,m— i,j—g,e 



dp, 



p 

nmgg 



i (Ar - X) {n - m) - -{n + m)j p^^^^ + Ka/ (n + l)(m + l)Pn+l,m+l,gg + 7l E Ppqeednpegdqmge 

p,q=0 

-^KT](^(n + m + 2 p^„,gg + a* Vn + lp^+i,™gg + a^Vm + lp^,„+i^gg) ^dt (A.5) 

+ (" + + l)Pn+U,n+l,gg + V" + IPn+l.mgg + OgVm + lPn,,n+l,gg + [l^gf - l) Pnrngg) (t) , 



dP«mee = -i (A^ - x) (?^ - m) - -(n + m) p^^^^ + k^J {n + l)(m + l)p^+i,™+i,ee + 1\P\ 



--\^ni^[{n^m + 2\ajf^ p^^^^ + a*\/n"+Tp^+i + aeV^^TMp^^^+i^ee) ^ dt (A.6) 

+ (y(n^TT)(^^rTT)p^+i^,„+i,ee + a*V?^+ lp^+i,™ee + + lp^^„+i^^p, + {^a^\^ - l) p^^^e) diV(0, 
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dp, 



nmeg 



(1 - t])k 



-i(^^a + {n- m)Ar + (n + m)x + ^ (ed/3*g + £*dPeg)^ ~ o'"'" (" + 



71 ^1 



7<^ - (1 - 



iIm(Q!*a 



3^)1 Pnmeg^^ 



+ (1 - ^?)'^ ( V + + l)Pn+l,m+l,e3 + /^eg + Ip^ ^^^^^^ - KgV^nTlPn + l,rneg) 



P 

nnieg 



(A.7) 



diV(t). 



I 

Noting that there are no mechanisms to populate To simplify this expression we look at the evolution of 
higher photon states in Eqs. (A. 6) and (A.7), and that the A's. Using that 
we have transformed to a frame where the initial pho- 
ton population is zero, we arrive at the equations for the 
diagonal elements of the qubit density matrix: 



dpf = dp^o,,, = ±7ipJJedi - \aif pfdt 



^pqeg 



+ [\a.^'-l)ptdN{t), (A.8) 

where the + (-) in the first term should be used when 
i ^ g (i ~ e). For the off-diagonal elements we have 



~PegV^d,p,q-i,egj 



^pqeg 



(A.ll) 



dt 



(Im(aga*)) = -- {edPlg + e*^l3eg) +2xRe(aga*) 



(A.9) 



-Kim {a*^ag) , 



(A.12) 



where 



^nmpqij = P^mij^pqij sxp (— ilm (aja*)) . (A. 10) and inserting the result from Eq. (A.8) gives us 



dA. 



nmpqeg 



i {uJa + B ~ r/Kim (alag)) — i{n — m)A,- — i(n + m)x 

(n- + "l) - y - 70 - + — |A>eg| j Kmpqegdt 



-(1 - T])k\/ (n + l)(m + l)Xn+l,m+l,pqegdt + {I - ri)KPegVm + lXn,m+l,pqegdt 
-(1 - ri)K(3*gVn+ lXn+l,mpqegdt + Peg y/pXnm,p-l,qeg dt 

n + m + \agf + |ae|^^ Kmpqeg + a^y/n + lXn+l,mpqeg + agVm + lXn,m+l,pqeg^ dt 



+ ^y(rr+T)(r7r+T)A„+i,m+i,pqeg -I- a*gVn + lA„+i^ 



mpqeg 



(A.13) 



Noting again that there are no mechanisms to populate trix: 

higher photon states, and that the initial photon popula- j_qb jt / -f- , tj r / * w Ti 

tion IS zero m the transformed frame, we get an equation " V ^ c j ^ y 2 

for the off-diagonal element of the reduced density ma- — r^; + — |/3 |^ ^p'^^dt 

-Y(l«.l' + Kl')p*di 

+ (aea;~l)ptdN{t). (A.14) 
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From Eqs. (A. 8), (A.14) and the identity 

= -77kP [n J (A. 15) 

we get the effective SME for photodetection 



dp' 



;qb 



B 



2 

70 + 



[cr^, p^'^] dt + 71 2? p'5'Mt(A.16) 



1 



where E[diV(i)] = 77K (nj;,!!^) dt. This is exactly what 
one gets by simply replacing a with IIq in all the mea- 
surement terms of the original SME. 

Normahzing Eq. (A.17) gives Eq. (13). Note that in 
the limit 77 — > Eq. (13) reduces to the effective ME 
without measurement. Secondly, averaging over an en- 
semble of quantum trajectories gives 



E 



-n>iM [n^n,] p'i'^dt + g [n„] p'5^d7V(t) 



r]KV [Ha] /o'^'Mt, 



(A.17) 



which lets us recover the effective ME without measure- 
ment. 
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